Mapping of quantum well eigenstates with semimagnetic probes 
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We present results of transmission measurements on CdTe quantum wells with thin semimagnetic 
Cdi-ajMn^Te probe layers embedded in various positions along the growth axis. The presence of the 
probes allow us to map the probability density functions by two independent methods: analyzing 
the exciton energy position and the exciton Zeeman splitting. We apply both approaches to map 
the first three quantum well eigenstates and we find that both of them yield equally accurate results. 
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I. INTRODUCTION 

The information about a quantum state is given by its 
energy E and its wave-function ip. Although it is easy to 
measure the former, the experimental access to the lat- 
ter is more difficult. In semiconductors, since the advent 
of epitaxial growth methods, band-gap engineering and 
tailoring of the eigenstates have become possible and re- 
trieving the information inscribed in the wave-functions 
became necessary. The access to -0 is gained through 
the probability density (PD) functions One way 

of probing of the PDs in quantum wells (QWs) is to in- 
troduce a highly localized potential perturbation with 
precisely controlled position along the growth axis. The 
perturbation shifts the eigenenergies of the system allow- 
ing to obtain the PD at the location of the probe in either 
an optical^ or a transport!^ experiment. Another way is 
to introduce a layer containing magnetic ions, which in 
an external magnetic field give rise to a Zeeman effect 
larger than in the case of an unperturbed QW— i^i^. This 
latter method allowed mapping of PDs in single^ or cou- 
pled multiple QWs^. In the first case, only the ground 
state PD was accessed and in the second, only relative, 
integrated PD values were obtained. In this work we 
apply both approaches to extract the PDs: we introduce 
magnetic probes and measure both the Zeeman effect and 
the shift of excitonic transitions in interband absorption. 
In this way, we map PD functions of the ground and the 
first two excited states of a CdTe QW sandwiched be- 
tween Cdi-yMgyTe barriers. As a PD probe, we use a 
layer of Cdi_a:Mna:Te, where a part of the Cd cations are 
substituted with Mn^+ ions. 

Incorporation of magnetic ions into a semiconductor 
matrix gives a new class of materials usually referred to 
as diluted magnetic semiconductors (DMS)^. Most com- 
monly, the substituting atoms are transition metal ions 
with partially filled d-shells (Mn^+ ions have a half filled 
d-shell), which gives rise to a locaHzed magnetic moment. 
Exchange interaction between localized spins of the d- 
shell electrons and band carriers leads to Zeeman effects 
enhanced by up to three orders of magnitude. To write 
the electronic wave-function in a DMS, we assume that 
the electrons adjust quasi-instantaneously to the arrange- 



ment of localized spins. In this adiabatic approximation 
the electronic wave-function reads^: 

*(r; Si, ...,Sn)^ *(r; S,) ^ ^{r; S^)<^>{{S^}) (1) 

where {Si} denotes the set of all quantum numbers de- 
scribing the system of magnetic ions. The s,p ~ d ex- 
change interaction is described by the Hamiltonian: 

Hsp-d = J'^'^ir- R,)S^<J (2) 
Ri 

where f and Ri are the spatial and a and Si are the spin 
coordinates of a band electron and a localized ion, re- 
spectively. As a consequence of the localized character 
of the c?-shell electrons, the exchange constant is usu- 
ally approximated by a collision term: J^'P~'^{f — Ri) = 

The s-d exchange interaction leads therefore to a con- 
duction band splitting given by: 

A£;, = S],($|5,|$)A^oa|0c(^»,>^^)n<^c(^.)l' (3) 

where [Xi,Yi,Zi), are the coordinates of a ith Mn ion, 
Nq is the number of cation sites per unit volume, and 
a is the s-d exchange integral. In the above, we fac- 
torized the electron wave-function into components de- 
pendent on the in-plane and perpendicular coordinates: 
-0(r) = (j)[x,y)(p{z). Such a procedure is not always 
justified when interband absorption is involved as the 
electron-hole Coulomb interaction mixes these degrees of 
freedom. However, we checked that in our case this mix- 
ing is negligible and in mapping experiments it leads to 
errors smaller than those resulting from compositional 
fluctuations and temperature instability. 

If the function ip{z) does not change substantially along 
the thickness of the probe layer, we can substitute the 
summation over Zi with a value of (p{z) at the Mn layer 
location Zmu- Moreover, assuming uniform distribution 
of Mn ions in the QW plane we can average 4'c{x, y) over 
in-plane ion coordinates and obtain the electron Zeeman 
splitting proportional to the layer magnetization M^—: 
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AE, = Noa\ipciZMn)\^\MX^,Y,)\^^,{^\S,\^) ^ 

^ Noa\^,{ZMn)? ■ Ml ^' 

where for Cdi_a;Mna:Te the conduction band exchange 
constant is^: Noa = 0.22 eV. 

It is thus seen from Eq. |4]that electron Zeeman spht- 
ting is proportional to the PD of finding an electron in the 
Mn layer. However, in an interband absorption experi- 
ment we measure the excitonic Zeeman splitting, which 
is a sum of electron and hole splittings: 

AEz = {Noa ■ \vc{ZMn)? - Nol3 ■ \ip4ZMn)\'') ■ Ml 

(5) 

where NqP = —0.88 eV is the valence band exchange 
constant^. 

It can be seen from the above that measuring excitonic 
Zeeman splitting for a series of samples, where the Mn 
ions are located at various positions Zmu , allows to map a 
PD function weighted with sp-d exchange integrals N^a 
and NqP contrary to the usual assumption^>^ that the 
heavy hole PD is mapped. 

Magnetic dopants not only give rise to magnetoopti- 
cal effects, but also introduce a local potential. In the 
first order perturbation theory, a potential of the form: 
V{S{z — Zmu)) located at the position Zmu shifts the 
electron energy of the eigenstate (pc{z) by: 

E',-E",^V\^c{ZMn)f (6) 

where V is the perturbing potential, given by the chemi- 
cal shift of the respective bands. Therefore, for electrons 
and holes the shift is proportional to conduction and va- 
lence band offset, respectively. As a result, the shift of 
the excitonic transition reads: 

AEs = E'x-E°x = gcV\^c{ZMn)\'' + gvV\(py{ZMn)f 

(7) 

where Qc and g-u are conduction and valence band offsets, 
respectively and E'^ is the energy of the unperturbed 
state. Therefore, measuring the shift of the exciton en- 
ergy, we can map a PD function weighted with band off- 
sets. In the following, we took a valence band offset = 
0.4i2iii and assumed a linear dependence of the chemical 
shift on the Mn composition^: V = 1592 meV • xmu- 



II. SAMPLES AND EXPERIMENT 

Designing samples for mapping experiments using 
semimagnetic probes one has to bear in mind that the 
profiles of CdTe/Cdi_i:Mna;Te interfaces are broadened 
along the growth axis due to complete exchange of Cd 
and Mn ions during growth and thus absence of segre- 
gation processesi^. Consequently, although we aim to 



obtain thin probe layers, the probe ions are always dis- 
tributed among a couple of adjacent monolayers with the 
composition profile peaked at Zmu- Moreover, we have to 
take into account the antiferromagnetic coupling between 
Mn ions, which decreases substantially the magnetization 
and as a result also the splitting, as seen in Eq. [H There- 
fore, the composition of the probe layers has to be low 
enough to assure a small number of nearest neighbor Mn 
pairs to avoid the antiferromagnetic coupling. 

The samples were grown on (001) oriented GaAs sub- 
strates by molecular beam epitaxy. Substrate temper- 
ature was 230° C, which assures a high sample qual- 
ity and relatively low interface broadeningi^. 3.5 ^m 
Cdi_yMgj,Te buffer was deposited before the growth of 
the QWs to relax the strain resulting from the lattice 
mismatch between the QW structure and the substrate. 
Next, five CdTe QWs, 117 A, wide were grown, sepa- 
rated by 300 A Cdi_j,Mgj,Te barrier layers. Magnesium 
composition y was chosen as high as 33% in order to 
assure that more than one confined state is present in 
the QW. Growth of each of the QWs was interrupted for 
the deposition of a single probe consisting of 2 monolay- 
ers of Cdi-i^Mn^^Te with intentional Mn molar fraction 
of 12%. For schematics of the probing heterostructure, 
see Fig. [Tb. Four samples with different positions of the 
probe layer along the QW axis were grown. Additionally, 
a reference sample with no probe was prepared. 

To obtain exciton transition energies and Zeeman split- 
ting we measured transmission as a function of magnetic 
field. To optically access the QWs, we first had to remove 
the nontransparent GaAs substrate, which was done by 
mechanical polishing and wet etching in hydrogen perox- 
ide. The thick transparent buffer produced Fabry-Perot 
oscillations, which obscured the absorption spectrum and 
thus most of it was removed by chemical etching in a 
0.6% solution of bromine in methanol. The sample was 
immersed in superfiuid liquid helium at a temperature of 
1.8 K. Magnetic field up to 7 T was appHed in Faraday 
configuration. A halogen lamp was used as a white-light 
source and the transmission signal was detected by a liq- 
uid nitrogen cooled CCD camera and a monochromator. 
In order to analyze transitions in two circular polariza- 
tions, a quarter wave-plate and a linear polarizer were 
placed in the way of the transmitted beam. 



III. RESULTS AND DISCUSSION 

Optical density spectra were evaluated according to 
Beer-Lambert law as ^ = — log(///o), where / and Iq 
are transmitted and incident beam intensities, respec- 
tively. In Fig. [H we show the spectra obtained from 
a sample with the probe layer located at the center of 
the QW. Transitions corresponding to three heavy hole 
excitons, labelled {nei,nhh) with n = 1,2,3 numbering 
electron and heavy hole states, can be resolved. Only 
diagonal transitions, i.e. those satisfying rigi = n^h, are 
observed. The oscillator strength of parity-allowed non- 
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diagonal transitions is very low, since in a deep rectan- 
gular QW the eigenstates are nearly orthogonal. The 
feature below the (2,2) transition is related to the H"^-" 
hole (1,1) exciton. The increase of the absorption ab 
2.1 eV is due to Cdi_yMgj,Te barrier excitons. Imi 
diately from Fig. 1 the effect of the Mn probe layer 
the exciton Zeeman splitting can be deduced: the O' 
number excitons exhibit a giant Zeeman effect since th 
is a non- vanishing PD of finding carriers in the cente] 
the QW (see Fig. 1). On the other hand, the Zeen 
splitting of the (2,2) state is limited to the direct inter 
tion between carriers and the magnetic field and so 
splitting is smaller then the transition line-width. 



meV. A very good agreement between the measured and 
calculated transition energies points out that the model 
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FIG. 1: Left: Optical density spectra for various magnetic 
fields obtained for a sample with the probe layer at the center 
of the quantum well. Right: Schematic of this sample shown 
together with the first three electron PD functions. 

The identification of the excitonic transitions in Fig. 
la is based on results of effective mass approximation 
calculations of electron and hole energy levels in the QW 
as a function of the magnetic field. In the calculations, 
we took into account the diffusion of the probe interfaces 
during growth. The broadened probe shape was modeled 
by a gaussian function^. The band-edge Zeeman split- 
ting of the probe was described by a modified Brillouin 
function^ with effective parameters and To refiecting 
the antiferromagnetic coupling between the Mn ions ad- 
justed to take into account the nonuniform number of 
nearest neighbors. Valence band states were calculated 
using a 4x4 Luttinger Hamiltoniani^. The lattice mis- 
match between the barrier and the QW layer introduced 
a strain, which was taken into account in the framework 
of the Bir-Pikus theory by adding a complete deforma- 
tion potential Hamiltoniaoi^. We neglected all the ex- 
citonic effects and used the exciton binding energy as a 
free parameter. In Fig. [21 we present the experimental 
and calculated transition energies for the same sample 
as in Fig. 1, i.e. with the Mn probe layer in the center 
of the QW. Heavy hole exciton binding energies result- 
ing from the presented fits were found to be between 14 
and 19 meV remaining in good agreement with calcula- 
tions in framework of an analytical model by Mathieu et 
ali^, which yields for (1,1) exciton a binding energy of 14 
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FIG. 2: Points: Exciton transitions measured for a sample 
with the probe layer at the center of the quantum well. Lines: 
The results of effective mass calculations allowing the identi- 
fication of the transitions. Full (empty) points correspond to 
transitions seen in (f ~) polarization. 

Using the above procedure for the reference sample 
with a fiat (i.e. without the probe layer) QW, we fit- 
ted the excitonic transitions and calculated electron and 
heavy hole PD functions |(/Jc(^)P and |(^„(z)|^. In Fig. 
[SK, we plot these functions weighted with exchange con- 
stants as derived in Eq. [5] for the first three QW eigen- 
states. On the same graph, we present the the exciton 
Zeeman splitting AEz measured as a function of the po- 
sition of the center of probe layer Zmu ■ A very good cor- 
relation between the Zeeman splitting and the PD value 
at Zmu is obtained confirming that the Zeeman split- 
ting provides a good estimation of the PDs. In Fig. [3b, 
we plot the fiat-QW PDs weighted with band offsets as 
derived in Eq. [7| On the same graph, we plot the ex- 
citon energy shifts AEs given by the difference between 
transition energies for the samples with probe layers and 
the reference sample with a fiat QW. The shifts for the 
(1,1) and (2,2) excitons are measured and for the remain- 
ing (3,3) is calculated since this transition could not be 
resolved in the reference sample. Again, a good correla- 
tion between the measured shift and the PD at the probe 
location is obtained. 

The accuracies of our mapping procedures are summa- 
rized in Figs. [3h and d, which show how the measured 
Zeeman splittings and excitonic shifts are correlated with 
values of PD functions weighted with exchange integrals 
and band offsets, respectively. In both cases, we obtain 
correlation coefficient values R « 0.9, proving a high ac- 
curacy of the approach. Inaccuracies are caused mainly 
by the fact that the probe layer has a nonzero thick- 
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FIG. 3: a) The PD functions of three lowest quantum well eigenstates weighted with exchange constants (lines) and corre- 
sponding Zeeman splittings (points, right scale) plotted as a function of the probe layer location, b) The PD functions of three 
lowest quantum well eigenstates weighted with conduction and valence band offsets (lines) and corresponding exciton shifts 
(points, right scale) plotted as a function of the probe layer location, c) and d) Correlations between weighted PD functions 
and Zeeman splitting and exciton shifts, respectively. 



ness, contrary to the assumption taken to derive Eqs. [5] 
and [71 Indeed, two monolayers deposited during growth 
are further broadened by the intermixing of the inter- 
face profile. To fit the excitonic transition dependencies 
on the magnetic field (see Fig. [2]), we had to introduce 
a gaussian broadened profiles with halfwidths between 
2.5 and 4 monolayers depending on the sample. Our ef- 
fective mass approximation calculations show that this 
nonzero thickness of the probe layers leads to a notice- 
able modification of the PD functions. In the PD func- 
tions compared with measured quantities in Figs. [3] also 
not included are the excitonic effects. The electron-hole 
Coulomb interaction modifies importantly the shape of 
PD functions with respect to a noninteracting case. A 
method to obtain Coulomb-correlated ^pc and ipv is based 
on a Hartree approach. In this calculation the electron 
wave function is self-consistently calculated by solving a 
one dimensional Schrodinger equation with an effective 
potential resulting from the hole wave function and vice- 
versaii^. Since in CdTe the hole is about 5 times heavier 
than the electron, the Coulomb-correlated Lpc tends to 
unperturbed ipv For this reason, previous works on Zee- 
man mapping successfully compared the results to only 
heavy hole PD functions^i^. 

Another issue not included in the above considerations 
is the variation of the exchange constants with confine- 
ment. As pointed out by Mackh et alt^ and Merkulov 
et alM, an admixture of higher fc- vector states leads to 
a significant decrease of the exchange parameters. Pos- 



sible reasons include turning on the kinetic exchange be- 
tween the conduction electrons and localized Mn ionsi^ 
and a hopping interference in the valence bancU^. In our 
case, both effects would have highest impact on the (3,3) 
state - the one with highest admixture of the nonzero-fc 
statesi^. We compared the PDs weighted with exchange 
integrals that were decreased according to the higher-fc 
states admixture and found a worse agreement with mea- 
sured Zeeman splittings. We therefore conclude that the 
reduction of the exchange integrals is less important than 
the effects related to a nonzero thickness of the probe lay- 
ers. 

In order to obtain a higher mapping accuracy one 
should design thinner (e.g. submonolayer) mapping 
probes that introduce a smaller perturbation potential. 
In that case, the observed Zeeman splitting and the shift 
of the excitonic transition will be smaller, but one will 
gain a weaker modification of the mapped PD functions. 
Another approach is to completely bury the semimag- 
netic probes in a CdMgTe QW exploiting the same band 
offsets for MgTe and MnTe with respect to CdTeifi>i^. 
In this method, the Mn composition of the probes and 
the Mg composition of the QW bottom are chosen to 
assure a flat QW potential so at zero magnetic fleld no 
perturbation is introduced. Applying a small magnetic 
field should make mapping feasible with only a minor 
modification of the PD functions. 

In summary, we have compared two independent meth- 
ods for mapping of the quantum well eigenstates with 
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semimagnetic probes. One is based on the analysis of 
the position of exciton transition which in the first or- 
der of perturbation theory is proportional to a value of 
a wave function weighted with band offsets. The sec- 
ond approach exploits the fact that for samples with thin 
semimagnetic probes the exciton Zeeman splitting is pro- 



portional to a value of a wave function weighted with ex- 
change constants. We find a good agreement between the 
calculated wave functions and measured excitonic posi- 
tions and Zeeman splittings and we conclude that both 
methods are equally well suited for mapping purposes. 
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